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This  paper  addresses  the  possible  connections  between  chaos,  the  unpredictable  behavior  of  solutions  of  finite  dimensional 
ssstems  of  ordinary  differential  and  difference  equations  and  turbulence,  the  unpredictable  behavior  of  solutions  of  partial 
differential  equations  It  is  dedicated  to  Martin  Kru.skal  on  the  occasion  of  his  60th  birthday 


I.  Introduction 

The  chaos  that  occurs  in  p.d.e.'s 
cannot  be  fathomed  by  legalesc 
so  we  apply  Occam's  razor 
and  using  a  laser 
study  structures  in  ring  Cannes 
An  appealing  idea  of  modern  dynamics  is  that 
the  complicated  and  apparently  stochastic  time 
behavior  of  large  and  even  infinite-dimensional 
nonlinear  systems  is  in  fact  a  manifestation  of  a 
deterministic  flow  on  a  low-dimensional  chaotic 
attractor.  If  the  system  is  indeed  low  dimensional, 
it  is  natural  to  ask  whether  one  can  identify  the 
physical  characteristics  such  as  the  spatial  struc¬ 
ture  of  those  few  active  modes  which  dominate  the 
dynamics.  Our  thesis  is  that  these  modes  arc  closely 
related  to  and  best  described  in  terms  of  asymptot¬ 
ically  robust,  mulnparamcter  solutions  of  the  non¬ 
linear  governing  equations.  We  find  it  hard  to 
define  this  robust  nature  precisely,  but  loosely 
speaking  the  idea  is  that  these  solutions  are  very 
stable  and  resilient  asymptotic  states.  They  may  be 
coherent  lumps  like  soluons  and  solitary  waves. 
They  may  have  the  form  of  coherent  wave  packets. 
They  may  have  self-similar  form.  They  need  not 
necessarily  be  the  asymptotic  states  which  develop 
as  t  tends  to  infinity;  structure-,  which  develop 

'Permanent  address  Dept  of  IMivmcv  Henoi-Uuit  Univer¬ 
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singularities  in  finite  time  like  those  involved  in  the 
collapse  of  Langmuir  waves  or  in  filamentation  in 
nonlinear  optics  are  also  candidates.  For  example, 
singular  solutions  of  the  Euler  equations  may  be 
useful  in  understanding  the  behavior  of  the 
Navicr-Stokes  equations  at  high  Reynolds  num¬ 
bers.  Singular  solutions  like  defects  and  disloca¬ 
tions  certainly  do  play  important  roles  in  the 
pattern  formations  arising  in  continuum  and  con¬ 
densed  matter  physics.  The  key  idea  is  that  each  of 
these  structures  is  a  natural  asymptotic  stale  that, 
by  virtue  of  the  various  force  balances  in  the 
governing  equations,  develops  an  identity  which 
does  not  easily  decay  or  disperse  away. 

One  can  envision  two  types  of  chaos  occurring. 
The  hrst  is  a  phase  or  weak  turbulence  which 
arises  when  there  is  an  endless  competition  be¬ 
tween  equally  resilient,  locahicd  coherent  struc¬ 
tures  which  are  infinite  time  asymptotic  states  and 
which  are  initiated  at  random  at  various  parts  of 
the  physical  domain.  Examples  of  this  type  of 
turbulence  are  solitary  waves  in  the  one-dimen¬ 
sional  complex  envelope  equation,  Raylcigh- 
Benard  roll  patterns  with  different  orientations 
and  the  oscillatory  skew  varicose  states  in  low 
Prandtl  number  convection.  It  is  to  be  expected 
that  such  dynamics  may  be  low  dimensional.  The 
second  type  of  chaos  is  much  more  dramatic  and, 
for  want  of  a  better  word,  may  be  described  as  an 
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Abstract 

A  new  mathematical  instability  technique  is  presented  and  applied  to 
determine  the  stability  properties  of  a  physically  important  class  of 
standing  waves  in  nonlinear  planar  optical  waveguides.  The  method  is 
illustrated  by  a  case  where  soliton  perturbation  techniques  or  variational 


methods  are  Inapplicable. 
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Extended  abstract 

Although  completely  integrable.  the  periodic 
sine-Gordon  (SG)  equation. 

+  sin  u  =  0. 

is  a  very  complicated  infinite  dimensional  Ham¬ 
iltonian  system.  However,  frequently  its  asymp¬ 
totic  states,  and  those  of  nearby  equations,  are 
dominated  by  only  a  few  degrees  of  freedom.  Here 
we  explain  the  origin  of  this  low  dimensionality 
and  give  a  recipe  for  determining  the  dominant 
degrees  of  freedom. 

We  consider  the  stability  of  an  V-phase  quasi- 
periodic  solution  u  to  (SG)  in  the  space  F  of  Cx. 
periodic  functions  of  x.  Linear  stability  is  treated 
in  [1.  2]  where  we 

(i)  construct  a  countably  infinite  collection  of 
solutions  of  the  linearized  problem  which  we 
conjecture  to  be  a  basis  for  the  function  space; 

(ii)  show  that  at  most  a  finite  number  of  elements 
(“modes")  in  this  collection  grow  exponen¬ 
tially  in  time  (the  rest  grow  ,u  worsi  poly  norm¬ 
ally): 

(iii)  provide  a  simple  criterion  whit.h  identities  the 
unstable  modes: 

(iv)  compute  the  growth  rates  .  ■!  these  unstable 
modes. 

This  work  (1|  shows  that  the  classical  modulational 


instability  is  merely  a  special  case  of  a  very  general 
instability.  Our  analysis  involves  spectral  theory 
and  complex  function  theory.  In  addition,  in  [1] 
we  use  geometric  methods  to  reveal  the  nonlinear 
phenomena  which  govern  the  origin  and  the 
saturation  of  the  instability.  We  think  that  this 
geometric  description  makes  the  qualitative  nature 
of  the  instability  completely  transparent. 

Our  analysis  proceeds  from  the  linear  spectral 
problem  which  can  be  used  to  integrate  (SG). 
Since  u(x.t)  is  periodic  in  x.  a  Floquet  spectral 
theory  is  appropriate.  The  cornerstone  of  this 
spectral  theory  is  the  Floquet  discriminant  J  :  F  * 
C  -*  C.  a  function  of  11  =  (u.  u,)e  F  and  £tC 
For  fixed  u  the  zeros  of  [J:(«.  £)  -  4]  are  the 
periodic  and  antipenodic  eigenvalues.  An  eigen¬ 
value  £(  is  called  a  simple  point  or  a  double  point 
according  to  its  order  as  a  zero  of  -  4.  Under 
the  sine-Gordon  How.  A(u.  E)  is  invariant  and 
contains  all  of  the  sine-Gordon  constants  of  mo¬ 
tion.  These' invariants,  together  with  a  collection  of 
dynamical  variables  ( pl }.  constitute  coordinates 
for  u.  The  dynamical  variables  {ji(}  are  in  1  1 
correspondence  with  the  critical  points  of  A.  and 
satisfy  O.D.E.'s  in  both  x  and  t. 

Our  analysis  of  this  instability  has  revealed 
several  features  [1]  of  the  geometry  of  isospectrul 
sets  which  were  previously  unknown  (at  least  to 
us).  Perhaps  the  most  striking  of  these  features  is 
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A  numerical  procedure  is  developed  which  measures  the  sine-Gordon  soliton  and  radiation  content  of  anv  field  (0.0.  i 
which  is  periodic  in  space  The  procedure  is  applied  to  the  held  generated  bv  a  damped,  dnven  sine-Gordon  equation  This 
field  can  be  either  temporally  periodic  (locked  to  the  dnven  or  chaotic.  In  either  case  the  numencal  measurement  shows  that 
the  spatial  structure  can  be  descnbed  by  only  a  few  spatiallv  localized  (soliton  wave-train)  modes  The  numencal  procedure 
quantitatively  identifies  the  presence,  number  and  properties  of  these  soliton  wave-trains.  For  example,  an  increase  of  span.il 
symmetry  is  accompanied  by  the  injection  of  additional  solitons  into  the  field 


1.  Introduction 

When  an  infinite-dimensional  dynamical  system 
(such  as  a  nonlinear  partial  differential  equation) 
is  strongly  perturbed,  a  rich  variety  of  distinct 
responses  is  observed.  These  include  temporal  and 
spatial  chaos,  temporal  chaos  with  spatial 
coherence,  and  intermittent  bursts  between  chaotic 
and  coherent  states.  We  are  particularly  interested 
here  in  coherent,  quasistable  spatial  structures  in  a 
chaotic  field  and  their  role  in  the  transitions  of  a 
field  into  and  out  of  chaos.  Physical  examples 
include  large  coherent  eddies  in  a  turbulent  flow, 
dislocations  in  Couette  flow,  and  solitary  profiles 
in  chaotic  waves. 

Support  is  gratefully  acknowledged  from  the  National  Sci¬ 
ence  Foundation,  the  U  S.  Department  of  Energy,  U  S  Air 
Force,  and  Los  Alamos  National  Laboratory 


The  latter  case  is  the  easiest  to  study.  It  can  be 
modeled  [1-4]  by  a  nonlinear,  dispersive  wave 
equation  in  one  spatial  dimension  which  is  per¬ 
turbed  by  dissipation  and  by  a  sinusoidal  driver. 
In  contrast  with  other  examples,  one  dimensional 
waves  (i)  are  relatively  easy  to  integrate  numeri¬ 
cally;  and  (ii)  possess  coherent  spatial  structures 
(solitary  waves)  which  are  very  well  understood  in 
the  absence  of  perturbations. 

A  concrete  example  is  the  driven,  damped  sine- 
Gordon  equation  with  periodic  boundary  condi¬ 
tions  of  period  L: 

<t>,,  ~  =  Tcosoit  -  a<t>,,  (1) 

e.4U -l,./)  =  e.*(«.ni  (l.hct 

<p(x.  t  =  0)  =  <#>,„(  x).  (lie  I 

<t>,(x.t  =  0)  =  uin(.x). 


0167-2789/86/S03.50  ©  Elsevier  Science  Publishers  B.V. 
(North-Holland  Physics  Publishing  Division) 


4 


SCT  FOR 
ERECTIONS 


Phystca  16D  (1986)  000-000 
North-Holland.  Amsitraim 


A  QUASI-PERiODIC  ROUTE  TO  CHAOS  IN  A  NEAR-INTEGRABLE  PDE 

J 

A.R.  BISHOP.  M.G.  FOREST*.  D.W.  McLAUGHLIN **  and  E.A.  OVERMAN  II* 

Theoretical  Division  ana  Center  for  Wonhnear  Studies.  Los  Alunios  Sationul  Laboratory.  Los  Alamos.  SM  8~ J95.  L  SA 


fce> 


Pattern  formation  and  iransiuons  to  chaos  are  described  (or  the  damped,  ac-dnven.  one-dimensional,  periodic  sine-Gordon 
equation.  In  a  nonlinear  Schrodineer  regime,  a  generic  quasi-penodic  route  to  intermittent  chaos  is  exhibited  in  detail  using  a 
range  of-dvnanucal  systems  diagnostics.  In  addition,  a  nonbnear  spectral  transform  is  exploited:  (i)  to  identify  and  quantify 
—  kcfepJitrtoaiiofbof  space-time  attractors  in  terms  of  a  small  number  of  sohton  modes'  of  the  underlying  integrable  systerfE  111 
to  use  these  anaivuc  coordinates  to  identify  ■■■^■■homoctinic  orbits  as  possible  sources  of  chat^^nd  mi)  to  demorNrfate 
the  sigmhcance  of  b»»«thee  wtimi  rnilnnmi  me  energy  transier  between  coherent  and  extended  states  in  this  cbaouc 
system. 


j 


1.  Introduction 

Solutions  of  nonlinear  evolution  equations  often 
exhibit  rich  patterns  in  space  and  time  which  can 
have  both  coherent  and  chaotic  components. 
Throughout  this  conference  it  has  been  apparent 
that  significant  progress  is  being  made  in  under¬ 
standing  such  ” space-time  complexity.”  We  have 
heard  about  parallel  developments  in  three  main 
areas:  (i)  partial  differential  equations  (1-3];  (ii) 
statistical  mechanics  [4];  and  (iiil  cellular  au¬ 
tomata  [5]  and  coupled  map  lattices  [6],  In  each  of 
these  three  areas  it  is  useful  to  distinguish  systems 
which  are  dissipation-dominated  from  nearly  con¬ 
servative  ones.  In  the  dissipation-dominated  case 
(e.g.  (3|),  only  a  few  stable  or  metastable  states 
(space-time  attractors)  exist,  and  the  organization 
into  patterns  occurs  because  the  dissipation  causes 
the  system  to  quickly  relax  to  these  states  The 
attractors  can  usually  be  identified  as  local  minima 
of  a  Lyapunov  free  energy  functional,  and  pattern 
-  selection  is  frequently  determined  from  considera¬ 
tion  of  boundary  conditions.  The  pnenomena  are 
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often  generic  for  relatively  large  classes  of  equa¬ 
tions  and  initial  data. 

In  near  conservative  cases,  the  organizational 
process  is  more  subtle.  Typically,  the  underlying 
conservative  problem  has  many  exact  solutions  x 
with  a  wide  variety  of  spatial  and  temporal  pat¬ 
terns.  Of  these,  some  will  resonate  with  the 


perturbation  and  adjust  to  stable  or  metastable 
solutions  of  the  near  conservative  system  (1.  2. 

7-9],  Thus,  the  organization  occurs  through  a 
nonlinear  resonance  process  [10. 11],  and  it  can  be 
quite  sensitive  to  both  the  particular  equation  and 
initial  data.  Nevertheless,  this  organization  can.  uu 
experimentally  substantial  parameter  regim^  j 
simplify  the  chaotic  responses  (^2r)  these  inrimtr  q  -f 
dimensional  systems,  (worn  il\ata  which  might  be 
anticipated  from  >iy  i  i  mu  i  iPi  I  I  i  '  mi  n 
sional  analogues  (e.g.  single  nonlinear  oscillators i  v. 

(7.  8].-Concemmg  chaos  in  coupled  oscillators.  \  .  .  C. 

more  is  twomlv  different!  -J — - — 

In  both  the  dissipation-dominated  and  near 

conservative  cases,  the  solutions  reside  in  an  in- - - 

finite  dimensional  phase  space  but  may  approach 
attractors  which  are  low  dimensional  (1.  7-9], 

Hence,  the  mathematical  techniques  developed  in 
modem  dynamical  systems  theory  [12]  are  hope¬ 
fully  relexant.  In  particular,  these  techniques 
explain  how  motion  near  a  low  dimensional  at- 
IM  PORT ANT 
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Abatract 

Ua  cooaldar  several  aaaaplaa  of  nonllnaar  avolutloo  aquatloaa  with  Initial 
data  that  ara  rapidly  oscillating  functions  of  tha  apaca  varlabla.  Wa  obtain  an 
affactlva  ayataa  of  nonllnaar  avolutlon  aquations  for  tha  various  momenta  of  tha 
aolutlon  by  a  aultlpla  aeala  aathod.  Wa  alao  show  how  In  ona  caaa  (tha  Carlanan 
nodal)  conpanaatad  compactness  gives  a  vary  general  way  of  obtaining  tha  affactlva 
aquatloaa  without  tha  uaa  of  aultlpla  scales. 

1  •  Introduction 

Wa  ara  lntaraatad  In  tha  bahavlor  of  aolutlona  of  dataralnlatlc,  nonllnaar 
avolutloo  aquatloaa  or  ayataaa  whan  tha  initial  data  ara  rapidly  oaclllatlng 
functions  of  tha  apaca  varlabla.  Such  problans  cannot  bn  solved  nuaarlcally  In  a 
dlract  way.  Thay  ara  fraquantly  analysed  by  obtaining  a  hierarchy  of  aquations  tor 
tha  various  nowante  of  tha  aolutlon.  This  hierarchy  la  Infinite  dua  to  the 
nonlinearity  of  tha  aquaclona  and  It  la  fraquantly  rendered  finite  by  various 
ad  hoc  cloaura  procedures. 

We  shall  cxaalne  closely  two  axaaplce  of  aquations  with  rapidly  oscillating 
data  In  tha  following  sections:  tha  Carlanan  aquations  and  tbs  groadwall  aquations. 
Thasa  exaaplas  ara  analysed  by  tha  usual  aayaptotlc  methods  of  aultlpla  scales  and 
aodulatlon  theory.  The  fora  of  tha  affective  equations  or  of  tha  hierarchy  la  than 
discussed  and  coapared  with  tha  hierarchy  obtained  directly.  In  tha  case  of  tha 
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Many-parameter  routes  to  optical  turbulence 
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The  output  of  an  externally  pumped  passive  nonlinear-optical  ring  resonator  is  known  to  exhibit 
chaotic  dynamics  under  variation  of  one  or  more  external  control  parameters.  The  dynamics  of  the 
internal  transverse  laser  beam  profile  is  described,  in  the  good-cavity  limit,  by  an  infinite¬ 
dimensional  discrete-time  map  in  function  space.  We  find  that  the  asymptotic  dynamical  behavior 
of  the  internal  resonator  field  is  in  marked  contrast  to  earlier  plane-wave  predictions.  Instead,  self- 
focusing,  self-defocusing  nonlinearities,  and  linear  diffraction,  combined  with  the  pump  and  feed¬ 
back  of  the  ring  resonator,  give  rise  to  strong  spatial  modulation  of  the  transverse  profile.  For  a 
self-focusing  nonlinearity,  spatially  coherent  robust  transverse  solitonlike  structures  persist  even 
though  the  beam  is  undergoing  temporally  chaotic  motion.  Using  attractor  embedding  techniques 
on  our  numerical  solutions,  we  isolate  distinct  routes  to  optical  turbulence  under  variation  of  a  num¬ 
ber  of  different  control  parameters.  In  particular,  we  provide  examples  of  Ruelle-Takens-Newhouse 
sequences,  intermittency,  and  a  nongeneric  bifurcation  involving  period  doubling  of  invariant  circles. 

A  main  result  of  the  paper  is  that  many  few-dimensional  attractors  can  coexist  in  different  regions 
of  the  infinite-dimensional  phase  space  at  fixed  values  of  the  external  control  parameters.  These  at¬ 
tractors  are  accessed  by  varying  the  systems  initial  conditions  and  can,  under  variation  of  an  ap¬ 
propriate  control  parameter,  independently  undergo  transition  to  chaos.  Interaction  between  neigh¬ 
boring  attractors  appears  to  be  responsible  for  numerically  observed  departures  from  generic 
behavior.  We  note  striking  similarities  between  some  of  our  numerically  generated  instability  se¬ 
quences  and  recent  experimental  observations  in  low-aspect-ratio  fluid  systems. 

MS  code  no.  AL3060  1986  PACS  numberis):  42.65.  -k,  42.60.Da,  42.3e.Tj 
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'  1  '  I  1  ! 

I.  INTRODUCTION  designing  low-aspect-ratio  experiments.  While  truncated 

models  can  exhibit  complex  bifurcation  structure,  the  na- 
The  transition  to  turbulence  in  physical  systems  is  an  ture  of  the  transition  to  chaos  becomes  sensitive  to  the 

old  subject  with  much  of  the  earlier  focus  being  on  fluid  level  of  truncation.11  There  is  in  fact  evidence  to  show 

dynamical  systems.  Understanding  fully  developed  tur-  that,  in  some  instances,  the  original  partial  differential 

bulence  and  its  onset  is  still  a  major  unsolved  problem,1  equations  from  which  these  models  are  derived  exhibit  no 

although  spectacular  success  has  been  achieved  over  the  turbulence  over  the  same  parameter  range, 

past  few  decades  in  understanding  the  onset  of  low-level  The  nonlinear-optical  model  considered  in  the  present 
turbulence  in  low-aspect-ratio  fluids.2  Experiments  on  paper  has  the  advantage  that,  unlike  the  Navier-Stokes 

Rayleigh- Bernard  convection3  and  Taylor-Couette  flow4  equations  of  fluid  dynamics,  the  full  partial  differential 

have  essentially  confirmed  the  mathematical  scenarios  equations  are  directly  amenable  to  numerical  study  over  a 

proposed  as  a  result  of  the  work  of  Newhouse,  Ruelle,  and  reasonably  wide  parameter  range.  Ad  hoc  truncations  to 

Takens,5  Feigenbaum,6  and  others.7  The  difficulty  with  plane  waves  or  TEM  resonator  modes  are  thereby  avoid- 

these  mathematical  models,  however,  is  that  they  rely  on  ed.  Our  results  show  that  such  truncations  are  inap- 

topological  ideas  or  on  the  study  of  the  universal  proper-  propriate  and  would  lead  to  erroneous  results.  The  nu- 

ties  of  few-dimensional  maps.  Bifurcation  parameters  ap-  merical  studies  show  that  low-level  turbulence  appears  to 

pearing  in  these  models  cannot  easily  be  identified  with  be  intrinsic  to  nonlinear-optical  systems,  unlike  fluids  as 

physical  stress  (control)  parameters  such  as  a  Reynolds  or  noted  above.  We  observe  a  close  analogy  between  our  nu- 

Prandtl  number  for  a  fluid  flow.  An  alternative  approach  merically  generated  optical  results  and  experimentally  ob¬ 
is  based  on  a  Galerkin  approximation,  whereby  the  served  instabilities  in  the  latter  low-aspect-ratio  fluid  sys- 

geometry  of  the  system  is  exploited  and  the  full  fluid  terns. 

equations  expanded  in  a  Fourier  senes.’  Truncation  of  Haken,10  in  1975,  was  the  first  to  point  out  the  close 
the  Fourier  series  to  a  few  spatial  modes  leads  to  a  system  similarity  between  the  Lorenz  equations  and  the  single- 

of  coupled  nonlinear  ordinary  differential  equations  that  mode  laser  equations  for  a  bad  cavity.  Originally,  the 

can  then  be  integrated  directly.  Such  a  procedure  led  physical  constraints  required  in  order  to  make  this  model 

Lorentz*  to  observe,  for  the  first  time,  chaotic  dynamics  valid  were  thought  to  be  too  severe  to  lead  to  a  successful 

in  a  coupled  system  of  three  ordinarv  differential  equa-  experimental  verification.  Recently,  however,  expen- 

tions.  However,  it  is  not  clear  to  what  extent  truncated  ments  using  far  and  mid  infrared  lasers11  suggest  that  the 

models  mimic  the  true  physical  situation  Cases  m  point  proper  bad-cavity  conditions  and  required  gain  can  be 

are  the  above  fluid  experiments  where  (he  number  of  de-  achieved  to  test  the  Haken  model.  No  direct  comparison 

grees  of  freedom  available  to  the  system  are  reduced  by  between  the  model  and  experiment  has  yet  been  made. 
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Numerical  studies  of  the  propagation  of  waves  guided  by  a  thin  film  bounded  by  two  self-focusing  media  reveal  a 
whole  new  class  of  nonstationary  solutions  tr.at  can.  in  their  own  right,  propagate  as  guided-wave  fields,  some  with 
solitonlike  components. 


Self-focusing  effects  in  slab-waveguiding  structures 
have  been  predicted  theoretically,1'5  and  in  a  few  cases 
observed  experimentally ,6"  to  lead  to  a  number  of 
novel  phenomena.  Guided-wave  dispersion  relations 
based  on  the  nonlinear  wave  equation  were  obtained 
initially  for  Kerr-type  media’*'5  and  subsequently  for 
more  general  nonlinearities.5  9  The  salient  theoretical 
predictions  were  that  both  the  propagation  wave  vec¬ 
tor  and  field  distributions  vary  with  guided-wave  pow¬ 
er  and  that  new  branches  with  power  thresholds  can 
exist.  These  in  turn  have  led  to  the  identification3*  of 
a  number  of  potential  all-optical  devices,  such  as  opti¬ 
cal  limiters,  lower-threshold  structures,  and  bistable 
and  switching  devices. 

The  stability  to  propagation  of  nonlinear  single- 
interface1112  and  thin-film  guided13-14  waves  has  been 
studied  recently.  For  nonlinear  waves  guided  by  a 
single  interface  between  a  self-focusing  medium  and  a 
power-independent  medium.  Akhmediev  er  a/. 1112 
found  the  solutions  to  be  stable  only  if  dP  c 3  >  0. 
where  P  is  the  guided-wave  power  and  J  is  the  effective 
index.  Although  similar  results  were  found  for  TE> 
waves  guided  by  thin  films  bounded  by  two  self-focus¬ 
ing  media,13-14  the  TE|  (Ref.  15)  case  was  more  com¬ 
plex.  In  this  Letter  we  report  further  numerical  in¬ 
vestigations  of  the  TE)  unstable  branches  that  reveal  a 
whole  new  class  of  nonstationary  guided  waves.  Un¬ 
der  appropriate  conditions,  these  solutions  exhibit  sol¬ 
itonlike  components  that  propagate  away  frc-m  the 
film  structure. 

This  study  deals  with  an  asymmetric  slab-wave- 
guide  geometry  for  which  all-optical  switching  and 
bistability  have  been  identified  as  possible  device  ap¬ 
plications.  ,y  The  TE0  field  is  written  as 

£(>,  n  =  z)exp[<(  -  „:j  -  c  c  .  1 1) 

with  k0  =  w/c;  3  is  the  effective  index.  3oth  cladding 
( -d  >  x)  and  substrate  (x  >  d)  layers  exhibit  a  positive 
(self-focusing)  nonlinearity  with  refractive  indices  giv¬ 
er.  by  n-TrJ  IV12)  =  n-\  +  «.]  H’!2,  y  =  c,s  (  r  >  d\  for 


the  cladding  and  the  substrate,  respectively,  with  n:  = 
ns  =  1.55  and  ac  -  2 a,.  The  film  medium  is  2d  thick, 
with  2ko d  =  12.2,  and  has  an  intensity-independent 
refractive  index  n/  =  1.57.  These  parameters  corre¬ 
spond  to  those  of  the  film  and  the  cladding  used  ex¬ 
perimentally  in  Ref.  6.  In  order  to  incorporate  propa¬ 
gation  effects,  we  must  study  a  mixed-type  linear- 
nonlinear  Schrodinger  equation,  which  is  more  com¬ 
plex  than  the  stationary  nonlinear  wave  equation.  It 
is  of  the  form 

2 i8k0  ^  =  -  ^  +  k02[32  -  n2(x.  |  W\ 2>)  U’  -  2 1 
dz  dx ' 

We  have  assumed  in  deriving  Eq.  (2)  that  the  field 
propagates  close  to  the  axis  of  the  guide  (along  c  i  and 
that  the  field  variation  along  the  z  axis  is  small  over 
one  wavelength.  Note  that  the  usual  stationary  non¬ 
linear  waves  are  z-independent  solutions  [H'(x,  cl  = 
IV’oU))  of  Eq.  (2).  The  properties  of  Eq.  (2)  are  well 
known  for  a  homogeneous  isotropic  nonlinear  Kerr 
medium,  and  a  whole  class  of  iV-soliton  and  radiation 
solutions  is  known  to  exist.16-1"  It  is  this  type  of  solu¬ 
tion  that  we  seek  in  the  present  case.  In  numerically 
solving  Eq.  (2)  we  have  used  the  beam-propagation 
method,1"  incorporating  a  split-step  fast  Fourier 
transform  with  512  transverse  sampling  points  and  a 
step  length  of  0.01  free-space  wavelength.  In  all  ex¬ 
amples  we  have  propagated  the  waves  a  minimum  or 
200  wavelengths  along  the  z  axis,  which  required  2-3  h 
of  CPU  time  on  a  Data  General  MV10000  computer 
Figure  1  summarizes  the  stability  properties  of  non 
linear  guided  waves  initially  launched  at  points  on  the 
nonlinear  dispersion  curves  obtained  for  the  guided 
wave  power  P  versus  effective  index  J.  Unstable  si  j 
tionary  waves  are  defined  as  waves  whose  field  d;s:r; 
butions  change  with  propagation  distance.  As  expect 
ed.  the  TEy-like  dispersion  curves  in  Fig.  1  show- 
similar  stability  properties  to  earlier  results1314  oh 
tained  for  2 kpd  =  1.256.  Loss  of  stability  always  oc¬ 
curs  on  the  upper  branch  (II),  except  for  the  region  ,n 
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The  stability  of  stationary,  TE„-type,  nonlinear,  thin-fiim  guided  waves  was  investigated 
numerically  for  both  symmetric  and  asymmetric  planar  waveguides  with  nonlinear  cladding  and 
substrate  layers.  It  is  found  that  large  regions  of  the  dispersion  curves  are  unstable  at  high  powers. 


Unique  properties'"7  have  been  predicted  for  waves 
guided  by  thin  films  when  one  or  more  of  the  guiding  media 
exhibit  a  field-dependent  refractive  index.  Self-focusing 
bounding  media  lead  to  multiple  new  branches  with  power 
thresholds,  as  well  as  field  distributions  whose  maxima  shift 
from  the  film  to  the  bounding  media  with  increasing  pow¬ 
er.1"7  In  fact,  this  geometry  has  been  identified  as  an  excel¬ 
lent  candidate  for  all-optical  switching,  with  our  without 
bistability  . "“7  To  date,  however,  theoretical  analysis  has  been 
based  solely  on  steady-state  solutions  to  a  nonlinear  wave 
equation  which  contains  an  intensity-dependent  refractive 
index.  The  salient  question  is  whether  these  wave  solutions 
are  stable  on  propagation,  and  the  consequences  of  possible 
unstable  regions  to  proposed  devices.  For  the  related  prob¬ 
lem  of  self-focusing  of  plane  waves  in  infinite  media.  Kolo- 
kolov’'  has  shown  that  the  solutions  are  stable  for  dP  /dfi  >  0 
where  Pis  the  power  and  [5  is  the  power-dependent  refractive 
index  Numerical  propagation  studies  of  nonlinear  waves 
guided  by  the  interface  between  a  self-focusing  and  a  power- 
independent  medium  by  Akhmediev  and  co-workers'1" 
have  led  to  a  similar  conclusion.  Recently,  a  theory  based  on 
phase  portraits  has  been  developed1 1  for  the  stability  of  TE„ 
waves  guided  by  a  thin  film  bounded  by  self-focusing  media, 
and  the  important  conclusion  of  that  work  is  that  the  waves 
are  unstable  on  negatively  sloped  branches  {dP/dP  < 0)  of 
the  nonlinear  dispersion  curve,  and  that  they  are  stable  on 
positively  sloped  regions  provided  that  self-focusing  occurs 
in  only  one  nonlinear  medium.  In  this  letter  we  report  a  test 
of  this  conclusion  via  a  numerical  investigation  of  the  stabil¬ 
ity  of  TE„  solutions  for  films  bounded  by  self-focusing  me¬ 
dia. 

The  geometry  analyzed  consists  of  a  film  ( jx|  <d.  re¬ 
fractive  index  n„ )  bounded  by  two  nonlinear  media  with  low 
power  indices  n ,  and  n:.  as  shown  in  Fig.  I  Since  the  numeri¬ 
cal  analysis  is  performed  in  the  slowly  varying  phase  and 
amplitude  approximation,  we  write  the  optical  field  as 

£(/-)  =  !.W'U.z +  c.c. ,  (1) 

where/?  is  the  effective  guided  wave  index.  k„  =  a >/c  and  the 
variation  in  the  amplitude  term  W(xj)  along  the  propaga¬ 
tion  direction  c  is  assumed  to  be  small  over  one  wavelength. 
The  refractive  index  in  the  various  media  is  given  by 


n:[\fV(x,z)\2]  =n 2  +ar\W(xj)\2 

T-  1.2  (|x|>rf) 

(2a) 

and 

n2[|  W'U,z)|J]  =  n2  (  |jc| <z/)  . 

(2b) 

Substituting  into  the  nonlinear  wave  equation  and  retaining 
only  the  first  derivative  of  fV{xj)  with  respect  to  z  leads  to 

dz  dx~ 

\{p2 -n2[\W(xj)\2]}W{xj)=0.  (3) 

Analytical  stability  analysis  of  Eq.  (3)  is  complicated  by 
the  fact  that  it  is  a  partial  differential  equation  and,  further¬ 
more,  is  a  Hamiltonian  system.  The  usual  stability  analysis 
for  dissipative  dynamical  systems  does  not  apply  ( unless  one 
deliberately  introduces  losses  into  the  problem).  One  reason 
for  the  difficulty  in  studying  the  stability  of  Eq.  ( 3 )  is  that 
many  of  the  eigenvalues  of  the  linearization  of  (3)  lie  on  the 
imaginary  axis  which  is  precisely  the  condition  for  instability- 
in  a  dissipative  dynamical  system.  Although  it  has  been  pos¬ 
sible1  1  to  perform  a  stability  analysis  forTEm  with  m  =  0,  it 
is  necessary  to  proceed  numerically  for  m^O.  First,  the 
steady-state  solutions  with  W'(xj)— W(x)  were  obtained  in 
the  usual  way  '"7  to  obtain  a  field  solution  corresponding  to  a 
particular  point  on  one  of  the  nonlinear  guided  wave  solu¬ 
tion  branches.  This  distribution  was  then  assumed  to  be 
launched  at  z  =  0  and  Eq.  (3)  was  solved  numerically  for 
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FIG.  1  Nonlinear  waveguide  geometry  studied  for  wave  stability. 
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Abstract:  Guided  waves  with  unique,  power-dependent  properties  arise 
when  one  or  more  of  the  media  bounding  a  guiding  film  exhibits  an 
intensity-dependent  refractive  index.  Previous  theoretical  work  on  this 
problem  has  been  based  formal  ism- 1 imited  to  Kerr-type  nonlinear  media  in 
which  the  change  in  refractive  index  is  quadratic  in  the  optical  field. 

In  this  paper  a  formalism  recently  reported  by  langbein  et  al  is  used  to 
investigate  nonlinear  guided  wave  solutions  in  more  realistic  material 
systems.  It  is  shown  numerically  that  saturation  of.  the  optically 
induced  change  in  the  refractive  index  can  dramatically  alter,  and  in 
some  cases  eliminate  the  more  interesting  power -dependent  features  of 
the  solutions.  Nonlinear  wave  solutions  are  also  investigated  for  a 
larger  class  of  media  that  are  characterized  by  refractive  indices  which 
depend  on  the  optical  field  raised  to  some  arbitrary  power. 
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The  nonlinear  waves  guided  by  the  interface  between  a  linear  and  a  non-Kerr-lilve.  nonlinear  medium  were  investigated 
numerically  using  the  formalism  recentlv  developed  bv  Langbein  et  al.  (Optics  Comm  5?  I19K5I  41*]  The  power  dependence 
of  the  propagation  constant  and  field  distributions  were  calculated  for  both  power-law  and  saturable  field-dependent  dielectric 
constants 


I .  Introduction 

It  has  been  shown  theoretically  that  the  interface 
between  two  dielectric  media  can  guide  TE  polarized 
waves,  provided  that  a  power  threshold  is  exceeded, 
and  that  one  of  the  two  media  exhibits  a  dielectric 
constant  which  depends  quadratically  (Kerr-like  me- 
diumlon  the  local  optical  field  [1-5],  To  date, 
these  investigations  have  been  based  on  the  nonlinear 
wave  equation  which  can  be  solved  analytically  [  1 ,6] 
for  a  relative  dielectric  constant  of  the  form  e  = 
ey  +  e^L(ZT)  with  =  oc^E2.  This  approach  has 

also  been  used  to  investigate  nonlinear  waves  guided 
by  thin  films.  Recently,  Langbein  et  al.  reported  {7 J 
a  new  technique  for  analysing  nonlinear  guided  wave 
phenomena  which  does  not  require  analytical  field 
solutions  to  the  nonlinear  wave  equation  in  order  to 
evaluate  the  power  dependence  of  the  propagation 
wavevector.  Their  approach,  which  requires  numerical 
integration  in  some  cases,  is  valid  for  arbitrarv  forms 
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of  eNL.  This  important  development  now  allows  real¬ 
istic  materials  to  be  considered  theoretically.  For  ex¬ 
ample,  all  materials  exhibit  a  saturation  in  the  opti¬ 
cally  induced  refractive  index  change  and  it  has  al¬ 
ready  been  shown  that  saturation  has  a  large  effect  on 
nonlinear  waves  guided  by  thin  dielectric  films  [8) . 
Furthermore,  the  dependence  of  the  dielectric  con¬ 
stant  on  the  field  is  not  quadratic  in,  for  example, 
semiconductor  materials  [9-1 1  ] .  In  this  paper  we  in¬ 
vestigate  the  nonlinear  waves  guided  by  a  single  inter¬ 
face  for  saturable  nonlinear  media,  and  non-quadratic 
power  law  media. 

2.  Review  of  general  theory 

The  details  of  the  technique  can  be  found  in  ref. 

[7 J  and  here  we  briefly  summarize  the  salient  points. 
Writing  the  TE  solution  fields  in  the  y'th  medium  in 
the  general  form 
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Abstract:  Selective,  efficient  excitation  of  TEo  nonlinear  guided  waves  is  demonstrated 

numerically.  For  a  thin  film  bounded  by  two  self-focusing  media,  three  different  field 
distributions  corresponding  to  the  same  power  can  be  excited  independently  by  suitable 


Gaussian  input  beams. 
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Abstract:  The  beam  propagation  method  is  employed  to  analyse  nonlinear  directional 
coupler  operation  for  various  combinations  of  nonlinear  materials  and  initially 
mismatched  guides.  A  mixed  focusmc  defocusing  configuration  is  found  to  give  optimal 
results. 
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Abstract:  We  demonstrate  numerically  that  external  excitation  of  a  nonlinear  waveguide 
can  produce  sequential  threshold  behavior  via  multi-soliton  emission  from  the  waveguide. 


This  behavior  is  similar  to  that  predicted  to  occur  at  a  nonlinear  interface. 
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Abstract:  Guided  waves  with  unique,  power-aependenT  properties  arise 
when  one  or  more  of  the  media  bounding  a  guiding  film  exhibits  an 
intensity-dependent  refractive  index.  Previous  theoretical  work  on  this 
problem  has  been  based  formal  ism- 1  imi  ted  to  Kerr-type  nonlinear  media  in 
which  the  change  in  refractive  index  is  quadratic  in  the  optical  fiela. 

In  this  paper  a  formalism  recently  reported  by  langbein  et  al  is  used  to 
investigate  nonlinear  guided  wave  solutions  in  more  realistic  material 
systems.  It  is  shown  numerically  that  saturation  of  the  optically 
induced  change  in  the  refractive  index  can  dramatically  alter,  and  in 
some  cases  eliminate  the  more  interesting  power-dependent  features  of 
the  solutions.  Nonlinear  wave  solutions  are  also  investigated  for  a 
larger  class  of  media  that  are  character ized  by  refractive  indices  wnich 
depend  on  the  optical  field  raised  to  some  arbitrary  power. 


